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Abstract. We consider the incremental computation of minimal unsat¬ 
isfiable cores (MUCs) of QBFs. To this end, we equipped our incremental 
QBF solver DepQBF with a novel API to allow for incremental solving 
based on clause groups. A clause group is a set of clauses which is in¬ 
crementally added to or removed from a previously solved QBF. Our 
implementation of the novel API is related to incremental SAT solving 
based on selector variables and assumptions. However, the API entirely 
hides selector variables and assumptions from the user, which facilitates 
the integration of DepQBF in other tools. We present implementation 
details and, for the first time, report on experiments related to the com¬ 
putation of MUCs of QBFs using DepQBF’s novel clause group API. 


1 Introduction 

Let ij) = Q. (j)he a QBF in prenex CNF (PCNF) where Q = QiXi... (5„x„ with 
Qi € {V, 3} is the prefix containing quantified propositional variables Xi and (j) 
is a quantifier-free CNF. Given a PCNF ip = Q.cj), an unsatisfiable core (UC) 
of '0 is an unsatisfiable PCNF 0' = Q' . such that Q' C Q and C 0. The 
prefix Q' is obtained from Q by deleting the quantified variables which do not 
occur in 0'. A minimal unsatisfiable core 0Mf7(0|^of 0 is an unsatisfiable core 
0' = Q'. 0' of 0 where, for every C G 0', the PCNF Q'. (0' \ {C}) is satisfiable. 

Incremental solving is crucial for the computation of MUCs in the context 
of propositional logic (SAT), e.g. |ll3l8ll3l23l24l28j . Modifications of a CNF by 
adding and deleting clauses in incremental solving are typically implemented by 
selector variables and assumptions 121911011 712012112212,51291 . An added clause 
C is augmented by a fresh selector variable s so that actually C U {s} is added. 
Via the solver API, the user assigns these variables as assumptions under which 
the CNF is solved to control whether a clause is effectively present in the CNF. 

* Supported by the Austrian Science Fund (FWF) under grant S11409-N23. We would 
like to thank Aina Niemetz and Mathias Preiner for helpful discussions. This article 
will appear in the proceedings of the 18th International Conference on Theory and 
Applications of Satisfiability Testing (SAT), LNCS, Springer, 2015. 

^ The terminology minimal unsatisfiable subformula (MUS) is equivalent to MUC. 







Different from the assumption-based approach, the SAT solver zChafQ [26] 
provides an API to modify the CNF by adding and removing groups (sets) of 
clauses. Clauses are associated with an integer ID of the group they belong to. 

In assumption-based incremental solving, clause groups may be emulated by 
augmenting all clauses in a group by the same selector variable. The user must 
specify the necessary assumptions via the API in all forthcoming solver invoca¬ 
tions to enable and disable the right groups. In contrast to that, zChaff allows to 
delete groups by a single API function call. In terms of usability, we argue that 
incremental solving by a clause group API is less error-prone, more accessible to 
inexperienced users, and facilitates the integration of the solver in other tools. 

We present a novel clause group API of our QBF solver DepQBF (version 4.0 
or laterin the style of zChaff. Different from zChaff, we implemented clause 
groups based on selector variables and assumptions to combine the conceptual 
simplicity of zChaff’s API with state of the art assumption-based incremental 
solving. As a novel feature of our API, the handling of selector variables and 
assumptions is entirely carried out by the solver and is hidden from the user. 
Our approach is applicable to any SAT or QBF solver supporting assumptions. 
Based on the novel clause group API of DepQBF, we implemented a tool to 
compute MUCs of PCNFs, a problem which has not been considered so far. 
Results on benchmarks used in the QBF Gallery 2014 illustrate the applicability 
of the clause group API for MUG computation of PGNFs. 


2 Implementing a Clause Group API 

DepQBF is a solver for PGNFs based on the QBF-specific variant of the DPLL al¬ 
gorithm [6] with learning [i2ii8i:nj . Since version 3.0 pnm] . DepQBF supports 
incremental QBF solving via an API to add and remove clauses in a stack-based 
way (cf. Fig. 3 in jUj). This API is suitable for solving incremental encodings 
where clauses added most recently tend to be removed again in subsequent solver 
calls, like reachability problems such as conformant planning m or bounded 
model checking mnj. The new clause group API of DepQBF, however, allows 
to add and delete clauses arbitrarily, which is necessary for the incremental com¬ 
putation of MUCs of PCNFs. We first present our novel approach to keeping 
selector variables invisible to the user, which is a unique feature of DepQBF. To 
this end, we distinguish between selector variables and variables in the encoding. 

Let S = (^ 1 ,... ,'0„) be a sequence of PCNFs. We consider variables over 
which the PCNFs ipi are defined as user variables because they are part of 
the problem encoding represented by S. When solving S incrementally, selector 
variables used to augment clauses in ipi are not part of the original encoding. 
Variables v are stored in an array VA indexed by an integer ID id{v) of v such 
that VA[id(v)] = v. User and selector variables reside in separate sections of VA: 

^ zChaff website (July 2015): https://www.princeton.edu/~chaff/zchaff.html 
^ DepQBF is free software: http://lonsing.github.io/depqbf/ 
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The total size of TT is ws. The user variable section has size us. The following 
invariants are maintained: VA[id(yy\ = v where id{v) < ms if r; is a user variable 
and us < id{v) < ms if m is a selector variable. If a new user variable v with 
id{v) > us is added via the solver API, then VA is resized together with the user 
variable section. In this case the selector variables are assigned new, larger IDs 
and copied to a new position in VA. Then the literals of selector variables are 
renamed according to the newly assigned IDs in all (learned) clauses and cubes 
present in the current PCNF in a single pass. Resizing only the selector variable 
section of VA does not require assigning new IDs to selector variables. Similar to 
implementations of other SAT or QBF solvers, the user is responsible to avoid 
unnecessarily large user variable indices and thus avoid resizing VA. 

The API of DepQBF prevents accessing selector variables in VA, which are 
hence invisible to the user. In contrast to traditional solver implementations, 
e.g. dni, where the user is responsible to maintain selector variables manually, 
the internal separation between user and selector variables allows to conveniently 
allocate and rename selector variables on the fly inside the solver and without 
any user interaction. This feature is particularly useful for solving dynamically 
generated sequences S = {tpi,.. . ,?/>„) of PCNFs where the exact user variable 
IDs in each tpi are unknown at the beginning. 

In the following, we present the novel clause group API of DepQBF along 
with the example shown in Fig. A new clause group is created by calling 
new_cls_grp() , which returns a unique unsigned integer cgid as the ID of the 
group. Each time a new group cgid is created, internally a fresh selector variable 
s is allocated in the array VA and associated with the group cgid. 

A group cgid must be opened by open_cls_grp(cgid) before clauses can be 
added to it. All clauses added via the API are associated with the currently 
opened group cgid by internally augmenting them with the selector variable s of 
group cgid. Groups must be closed by close_cls_grp(cgid) . When solving the 
current PCNF by sat (), internally the selector variables of all created groups are 
assigned false as assumptions to effectively activate the clauses in these groups. 

Deleting a group by delete_cls_grp(cgid) invalidates its ID. When solving 
the current PCNF by sat (), internally the selector variables of all deleted groups 
are assigned true as assumptions to deactivate the clauses in all deleted groups 
and all learned clauses derived therefrom. Deleted clauses are physically removed 
from the data structures in a garbage collection phase if their number exceeds 
a certain threshold. Clauses which are added to the PCNF without opening a 
group by open_cls_grp(cgid) before are permanent and cannot be deleted. 

In contrast to deletion, clause groups can also be deactivated by calling 
deactivate_cls_grp(cgid) . When solving the current PCNF by sat(), in¬ 
ternally the selector variables of deactivated groups are assigned true simi¬ 
larly to deleted groups. However, clauses in deactivated groups are never re¬ 
moved from the data structures. Deactivated groups are activated again by 

















int main (int argc, char ** argv) {. 


...//continued from left column. 
Result res = sat(s); 
assert(res == RESULT_UNSAT); 
ClauseGroupID *rgrps = 


Solver *s = createO; 


new_scope_at_nesting 
(s,qTYPE_FORALL,l); 


add(s, 1 );add(s, 2 );add(s, 0 ); 


get_relevant_cls_grps(s); 
assert(rgrps [ 0 ] == id 2 ); 
reset(s); 


new_scope_at_nesting 
(s,qTYPE_EXISTS,2); 


add(s,3);add(s,4);add(s,0); 


ClauseGroupID idl = new_cls_grp(s); 
open_cls_grp(s,idl); 
add(s,-l);add(s,-3);add(s,0); 
close_cls_grp(s,idl); 


deactivate_cls_grp(s,rgrps[ 0 ]); 
res = sat(s); 

assert(res == RESULT_SAT); 
reset (s); 


ClauseGroupID id2 = new_cls_grp(s); 
open_cls_grp(s,id2); 
add(s,1);add(s,2);add(s,4);add(s,0); 
add(s,1);add(s,-4);add(s,0); 
close_cls_grp(s,id2); 

...//continues on right column. 


activate_cls_grp(s,rgrps[ 0 ] ); 
free(rgrps); 


delete_cls_grp(s,idl); 
res = sat(s); 

assert(res == RESULT_UNSAT); 
delete(s); } 


Fig. 1. Clause group code example. Variables Xi are encoded as integers i. Given the 
PCNF -i/i := yxi,X 2 Jx 3 , X 4 . C'iAC' 2 AC '3 with Ci = (-'ri V-'Xs), C 2 = (X1VX2VX4), C3 = 
(xi V “’*4), Gi is put in group idl and C2,C3 in group id2. An unsatishable core 
consisting only of group id 2 is extracted from ip. Deactivating group id 2 results in the 
PCNF VsiBxs. Cl. Activating id 2 again and deleting idl yields Vri, X 2 ^X 4 . G 2 A G 3 . 


activate_cls_grp(cgid). Selector variables of activated groups are assigned 
false when solving the current PCNF. 

DepQBF also allows for traditional incremental solving where the user han¬ 
dles selector variables manually nnj. Implementations of this approach like Min- 
iSAT, for example, allow to physically delete clauses by first adding a unit clause 
containing a selector variable and then simplifying the formula based on unit 
clauses. This is in contrast to DepQBF where the formula is not simplified based 
on unit clauses to avoid the internal elimination of variables, which may be 
unexpected by the user. 

If the current PCNF has been found unsatishable by sat(), then calling 
get_relevant_cls_grps 0 returns an array of the IDs of those groups which 
contain clauses used by the solver to determine unsatishability. The clauses in 
these groups amount to an unsatishable core of the PCNF. That core is obtained 
by internally collecting all selector variables relevant for unsatishabilitjj^ and 
mapping them to the respective clause group IDs. 


Similar to the function analyzeFinal in Minisat, for example. 
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QBF Gallery Track 

#m 

ut 

mt 

\CNF\ 

\MUC\ 

#c 

r r 

applications (190 of 735): 

182 

6,304 

7,941 

4,744,494 

73,206 

81,631 

6.1% 2.9% 

QBFLIB (58 of 276): 

46 

1,009 

2,264 

323,497 

34,777 

36,888 14.1% 5.1% 

preprocessing (38 of 243): 

34 

1,623 

1,080 

451,197 

23,220 

24,572 

4.0% 2.2% 


Table 1. Statistics for unsatisfiable instances from the QBF Gallery 2014 where MUCs 
were successfully computed. Numbers of solved instances out of total ones are shown 
in parentheses. MUCs computed {#m), total time to solve the initial unsatisfiable 
instances (ut) and to compute the MUCs (mt), total number of clauses in initial for¬ 
mulas (I CNF\) and in MUCs {\MUC\), total number of QBF solver calls (#c), and the 
average (r) and median (f) sizes of MUCs relative to the respective CNF sizes. 


3 Computing Minimal Unsatisfiable Cores of QBFs 

In contrast to theory [16], the computation of MUCs of PCNFs in practice has 
not been considered so far. Approaches to nonminimal UCs of PCNFs were 
presented in the context of checking Q-resolution refutations of PCNFs (SO] and 
QMaxSAT [U] . For the first time we report on experiments related to the compu¬ 
tation of MUCs of PCNFs. To this end, we implemented a tool to incrementally 
compute MUCs of PCNFs using the clause group API of DepQBF as follows. 

Given an unsatisfiable PCNF ipQ = Q. cj), first every single clause of t/jQ is 
put in an individual clause group. Let '0 := i/jq. The PCNF ip is solved and 
a UC Ip' = Q' . (p' is extracted by get_relevant_cls_grps. Then ip is replaced 
by Ip' by deleting the clause groups which do not belong to ip' from ip. Given 
the updated ip = Q.cp, every clause C G (p is checked by solving the PCNF 
Ip" = Q. ((/> \ {C}). To this end, the group containing C is deactivated. If ip" 
is satisfiable then C is part of an MUC and hence C is activated again iC is a 
transition clause jUj). Otherwise, a UC ip' of ip" is extracted, ip is replaced by 
the UC Ip' like above, and again every clause in the updated ip is checked. After 
every clause in the current ip has been checked, the final ip is an MUC of ipQ. The 
number of solver calls in this well-known elimination-based algorithm is linear in 
the size of ipo 11.112,1I24| . It applies iterative clause set refinement |^l51?r| by UCs. 
UCs are extracted by selector variables [T] in get_relevant_cls_grps, which is 
in contrast to extraction based on resolution proofs [27128]. The algorithm is 
common to compute MUCs of CNFs but has not been applied to PCNFs so far. 

Using our tool, we computed MUCs of instances from the applications (AT), 
QBFLIB (QT), and preprocessing (PT) tracks of the QBF Gallery 2014^ We 
preprocessed the instances from AT and QT using Bloqqer [5]. In total, we 
allowed 900s of wall clock time and seven GB of memory to solve an instance 
by DepQBF and to compute an MUC. Table [^ summarizes the results of our 
experiments run on an AMD Opteron 6238 at 2.6 GHz under 64-bit Linux. 
MUCs were successfully computed for 95% of the solved unsatisfiable instances 
in AT (79% of QT and 89% of PT). On average, MUC computation took 43s 


® http://qbf.satisfiability.org/gallery/ 














in AT (49s in QT and 31s in PT). When increasing the total timeout to 3600s, 
then 186 MUCs were computed in AT (48 in QT and 36 in PT). 

Iterative clause set refinement by UCs potentially reduces the number of 
solver calls. In the worst case, there is one solver call per each single clause in 
the initial PCNF '0o- However, on average there was one solver call per 58, 8, 
and 18 clauses in AT, QT, and PT, respectively. 

The physical deletion of clauses not belonging to a MUC reduces the mem¬ 
ory footprint and the run time. The plot below shows the sorted total run 
times (y-axis) of the MUC workflow on instances in TT where MUCs were suc¬ 
cessfully computed (x-axis). If clauses are deleted by delete_cls_grp (UC-d) 
then 182 MUCs are computed but only 169 if clauses are permanently deac¬ 
tivated by deactivate_cls_grp instead (UC-nd). We attribute this effect to 
overhead caused by deactivated clauses 
still present in the data structures. Only 
79 MUCs are computed without itera¬ 
tive clause set refinement by UCs us¬ 
ing get_relevant_cls_grps and instead 
checking every clause in ^j;Q one by one 
(OBO). We made similar observations for 
QT and PT. On instances where an MUC 
was computed by both UC-d and UC-nd, 
in general UC-nd is slower (up to -1-316% 
on PT) and has a larger memory footprint 
(up to -1-70% on AT). The difference between UC-d and OBO is more pro¬ 
nounced, where in general OBO is slower (up to -1-4126% on PT) and has a 
larger memory footprint (up to -1-243% on AT). 

Our experiments show that physical deletion of clauses by delete_cls_grp 
(UC-d) and the extraction of UCs by get_relevant_cls_grps based on selector 
variables are crucial for the computation of MUCs of PCNFs. These features are 
provided directly by the novel clause group API of DepQBF. 
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4 Conclusion 

We presented a novel API of our solver DepQBF for incremental QBF solving 
based on clause groups and its application to MUC computation. The clause 
group API is conceptually simple yet employs state of the art approaches to 
assumption-based incremental SAT solving. Improvements of assumption-based 
incremental solving |2I17I29| are also applicable to our implementation. 

The API encapsulates the handling of selector variables and assumptions 
entirely inside the solver. This is a unique feature of DepQBF, which facilitates 
its integration in other tools. It is particularly useful for solving dynamically 
generated sequences of PCNFs where the exact variable IDs are unknown at the 
beginning. The clause group API is general and fits any search-based SAT and 
QBF solver capable of solving under assumptions. 





A potential application of the clause group API is (M)UC extraction of PC- 
NFs in core-guided QMaxSAT [2] and SMT, similar to SAT-based UC extrac¬ 
tion in SMT [7]. Further, our API readily supports the extraction of high-level 
UCs [19127128] where, different from our experiments with MUC computation, 
multiple clauses are put in a clause group. We applied the novel clause group 
API of DepQBF to compute MUCs of PCNFs for the first time. Our results 
indicate the efficiency and applicability of our implementation. As future work, 
we want to integrate incremental preprocessing in DepQBF in a way where the 
implementation details are hidden by the API [5S]. 
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A Additional Experimental Data 


QBF Callery Track 

#m 

ut 

nit 

\CNF\ 

|MUC| 

#c 

r f 

applications (190 of 735): 

169 

4448 

7346 

1,604,908 

46,025 

53,080 

6.4% 3.0% 

QBFLIB (58 of 276): 

46 

982 

3,308 

323,497 

34,832 

37,020 14.1% 5.2% 

preprocessing (38 of 243): 

31 

613 

1658 

368,002 

14,673 

15,913 

3.4% 2.2% 


Table 2. Related to Table[^(same column labels): computation of MUCs where groups 

of clauses not being part of a MUC are deactivated instead of deleted (“UC-nd” in the 
plot in Section]^. Except for the QBFLIB track, fewer MUCs were computed when 
deactivating clauses than when deleting clauses. 


QBF Gallery Track 

#m 

ut 

mt 

\CNF\ 

\MUC\ 

#c r f 

applications (190 of 735): 

79 

238 

12,004 

310,986 

22,367 

311,065 10.6% 3.8% 

QBFLIB (58 of 276): 

34 

608 

5,970 

169,853 

20,463 

169,887 17.4% 7.9% 

preprocessing (38 of 243): 

21 

270 

5,326 

175,162 

5,233 

175,183 3.2% 2.6% 


Table 3. Related to Table (same column labels): computation of MUCs without 
iterative clause set rehnement by UCs but instead checking every clause in the initial 
PCNF one by one (“OBO” in the plot in Section]^. This approach performs consid¬ 
erably worse than either variant of iterative clause set refinement by UCs (UC-d and 
UC-nd in Tables and [^ . 



#m 

mt 

mem 

\CNF\ 

\MUC\ 

#c 

r 

r 

UC-d: 

46 

2,264 

2,175 

323,497 

34,777 

36,888 

14.1% 

5.1% 

UC-nd: 

46 

3,308 

2,948 

323,497 

34,832 

37,020 

14.1% 

5.2% 


Table 4. Related to Table (same column labels except mem which is the total 
amount of memory) and to the plot in Section]^ comparison of the MUC workflow 
on instances from the QBFLIB track where a MUC was computed by both iterative 
clause set refinement by UCs where clauses are deleted (UC-d) and deactivated (UC- 
nd), respectively. Whereas the sizes of MUCs and numbers of solver calls are similar for 
UC-d and UC-nd, UC-nd is slower (-1-46%) and has a larger memory footprint (-1-35%). 















#m 

mt 

mem 

\CNF\ 

\MUC\ 

#c 

r 

r 

UC-d: 

31 

398 

1,672 

368,002 

14,586 

15,785 

3.4% 

2.2% 

UC-nd: 

31 

1,658 

1,988 

368,002 

14,673 

15,913 

3.4% 

2.2% 


Table 5. Like Table but for instances from the preprocessing track. UC-nd is slower 
(+316%) and has a larger memory footprint (+18%). 


#m 

mt 

mem 

\CNF\ 

\MUC\ 

#c 

r 

r 

UC-d: 169 

4,118 

5,155 

1,604,908 

44,979 

52,049 

6.4% 

2.9% 

UC-nd: 169 

7,346 

8,779 

1,604,908 

46,025 

53,080 

6.4% 

3.0% 


Table 6. Like Table but for instances from the applications track. UC-nd is slower 
(+78%) and has a larger memory footprint (+70%). 



#m 

mt 

mem 

\CNF\ 

\MUC\ 

#c 

r 

f 

UC-d: 

34 

1,235 

1,040 

169,853 

19,603 

21,192 

15.7% 

5.1% 

OBO: 

34 

5,970 

1,743 

169,853 

20,463 

169,887 

17.4% 

7.9% 


Table 7. Related to Table (same column labels except mem which is the total 
amount of memory) and to the plot in Section comparison of the MUC workflow 
on instances from the QBFLIB track where a MUC was computed by both iterative 
clause set refinement by UCs where clauses are deleted (UC-d) and without iterative 
clause set refinement and checking clauses one by one instead {OBO), respectively. In 
addition to the solver calls for solving the initial unsatishable PCNFs, with OBO there 
is one solver call for each clause in the initial PCNF. OBO is slower (+383%) and has 
a larger memory footprint (+67%). 


#m mt 

mem 

\CNF\ 

\MUC\ 

#c 

r 

r 

UC-d: 21 126 

850 

175,162 

4,977 

5,861 

3.2% 

2.2% 

OBO: 21 5,326 

1,135 

175,162 

5,233 

175,183 

3.2% 

2.6% 


Table 8. Like Tablebut for instances from the preprocessing track. OBO is slower 
(+4126%) and has a larger memory footprint (+33%). 


#m 

mt 

mem 

\CNF\ 

\MUC\ 

#c 

r 

r 

UC-d: 79 

1,724 

839 

310,986 

19,878 

23,665 

10.3% 

2.9% 

OBO: 79 

12,004 

2,879 

310,986 

22,367 

311,065 

10.6% 

3.8% 


Table 9. Like Table [7| but for instances from the applications track. OBO is slower 
(+596%) and has a larger memory footprint (+243%). 
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Fig. 2. Like the plot in Sectionbut for instances from the QBFLIB track. 
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Fig. 3. Like the plot in Section]^ but for instances from the preprocessing track. 







